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IMAGE DATA COMPRESSION EMPLOYING MULTIPLE 
COMPRESSION CODE TABLES 



FIELD OF THE INVENTION 

5 The present invention relates generally to a field of image compression and 

decompression, and more particularly to a technique for rapidly and optimally 
compressing and decompressing image data through the use of one or more 
compression code tables selected from a family of predefined tables 



1 0 BACKGROUND OF THE INVENTION 

Various techniques have been proposed and are currently in use for 
analyzing and compressing large data files, such as image data files. Image data 
files typically include streams of data descriptive of image characteristics, typically 
of intensities or other characteristics of individual pixel elements or pixels in a 

15 reconstructed image. In the medical field, for example, large image files are 

typically created during an image acquisition or encoding sequence, such as in an x- 
ray system, a magnetic resonance imaging system, a computed tomography imaging 
system, and so forth. The image data is then processed, such as to adjust dynamic 
ranges, enhance certain features shown in the image, and so forth, for storage, 

20 transmittal and display. 



While image files may be stored in raw and processed formats, many image 
files are quite large, and would occupy considerable memory space. The increasing 
complexity of imaging systems also has led to the creation of very large image files, 
25 typically including more data as a result of the usefiil dynamic range of the imagmg 

system and the size of the matrix of image pixels. 



hi addition to occupying large segments of available memory, large image 
files can be difficult or time consuming to transmit from one location to another, hi 
30 a typical medical imaging application, for example, a scanner or other imaging 

device will typically create raw data which may be at least partially processed at the 
scanner. The data is then transmitted to other image processing circuitry, typically 
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including a programmed computer, where the image data is further processed and 
enhanced. Ultimately, the image data is stored either locally at the system, or in a 
picture archiving and communications system (PACS) for later retrieval and 
analysis. In all of these data transmission steps, the large image data file must be 
5 accessed and transmitted from one device to another. 

Compression techniques have been developed that apply various algorithms 
and approaches to conversion of original image data to a compressed form for 
transmission and storage. One such approach is based upon assignment of 

10 compressed data code by reference to a table, commonly referred to as a 
compression table. This approach is based on the probability (or the fi-equency) of 
occurrence of different levels, typically gray levels or intensity levels, for each pixel 
in an image, represented by corresponding binary values in the image data stream, 
hi general, compression code table permits more fi-equently occurring values to be 

15 assigned a shorter compressed data code than less frequently occurring values. 

Compression ratios in such techniques may, how^ever, be highly dependent upon the 
relative frequencies of occurrence of the values across the dynamic range of the 
image data. 

20 There is a present need for an improved technique for compressing and 

decompressing image data which provided higher relative compression ratios in a 
computationally efficient manner. There is a particular need for a technique which 
can be apphed to new and existing compression systems, and which can be adapted 
to various systems depending upon the characteristics of the images to be handled. 

25 

SUMMARY OF THE INVENTION 

The present invention provides an image data compression and 
decompression technique designed to respond to these needs. The technique is 
applicable to a wide variety of imaging fields, and is particularly well suited to 
30 medical diagnostic imaging systems. Compression in accordance with the technique 

is based upon reference to a set of compression tables, such as compression code 
tables which are predefined to accommodate different levels to data variation or 
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entropy in the image data stream. The compression code tables are stored for 
reference during the compression process, and compression performance is 
evaluated by application of the candidate tables to determine which table or tables 
provide the best compression. 

5 

In a presently preferred embodiment, the family of compression code tables 
offer differing compressed data code and differing compression ratios. The image 
data stream is analyzed in subregions to identify which compression code table 
provides the best compression for each subregion. The appropriate code table is 

1 0 selected in accordance with this analysis, and several such tables may be selected for 

different regions of the image data. Higher entropy regions, for example, may 
require application of tables providing somewhat reduced compression ratios, while 
lower entropy regions may use compression code tables providing much higher 
compression rations. The resulting overall compression ratio is therefore improved 

15 as compared to techniques employing a single table for an entire image. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a diagrammatical representation of a picture archiving and 
communication system or PACS for receivuig and storing image data in 
20 accordance with certain aspects of the present technique; 

Fig. 2 is a diagrammatical representation of contents of a database for 
referencing stored image data in files containing multiple image data sets, 
compressed data, and descriptive information; 

Fig. 3 is a representation of a typical image of the type received, 
25 compressed, and stored on the system of Fig. 1 ; 

Fig. 4 is a graphical representation of intensities of pixelated data across an 
image, subdivided into subregions for compression of the subregions optimally 
based upon characteristics of the subregions; 

Fig. 5 is a diagrammatical representation of a pixel neighborhood used in 
30 the analysis of the image data for compression purposes; 
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Fig. 6 is a flow chart illustrating exemplary control logic for compressing 
and decompressing image data in accordance with aspects of the present 
technique; 

Figs. 7, 8, 9, 10, 11 and 12 are reference lookup tables in the form of 
5 compression code tables used to optimally compress subregions of image data in 

accordance with the present technique during the process illustrated in Fig. 6; and 

Fig. 13 is a diagrammatical representation of an exemplary image data set, 
including a descriptive header, a compression header, and blocks of compressed 
data by subregion. 

10 

DETAILED DESCRIPTION OF THE INVENTION 

Fig. 1 illustrates a picture archive and communication system or PACS 10 
for receivmg, compressmg and decompressing image data. In the illustrated 

15 embodiment, PACS 10 receives image data from several separate imaging systems 

designated by reference numerals 12, 14 and 16. As will be appreciated by those 
skilled in the art, the imaging systems may be of various type and modahty, such as 
magnetic resonance imaging (MRI) systems, computed tomography (CT) systems, 
positron emission tomography (PET) systems, radio fluoroscopy (RF), computed 

20 radiography (CR), ultrasound systems, and so forth. Moreover, the systems may 
include processing stations or digitizing stations, such as equipment designed to 
provide digitized image data based upon existing fihn or hard copy images. It 
should also be noted that the systems supplying the image data to the PACS may be 
located locally with respect to the PACS, such as in the same institation or facility, 

25 or may be entirely remote from the PACS, such as in an outlying clinic or affiliated 

institution. In the latter case, the image data may be transmitted via any suitable 
network link, including open networks, proprietary networks, virtual private 
networks, and so forth. 

30 PACS 10 includes one or more file servers 18 designed to receive and 

process image data, and to make the image data available for decompression and 
review. Server 18 receives the image data through an input/output interface 19. 
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Image data may be compressed in routines accessed through a 
compression/decompression interface 20. As described more fully below, interface 
20 serves to compress the incoming image data rapidly and optimally, while 
maintaining descriptive image data available for reference by server 18 and other 
5 components of the PACS. Where desired, interface 20 may also serve to 

decompress image data accessed through the server. The server is also coupled to 
internal clients, as indicated at reference numeral 22, each client typically including 
a work station at which a radiologist, physician, or clinician may access image data 
fi-om the server, decompress the image data, and view or output the image data as 
10 desired. Chents 22 may also input information, such as dictation of a radiologist 

following review of examination sequences. Sunilarly, server 18 may be coupled to 
one or more interfaces, such as a printer interface 24 designed to access and 
decompress image data, and to output hard copy images via a printer 26 or other 
peripheral. 

15 

Server 28 may associate image data, and other work flow information within 
the PACS by reference to one or more file servers 18. In the presently contemplated 
embodiment, database server 28 may include cross-referenced information regarding 
specific image sequences, referring or diagnosing physician information, patient 

20 information, background information, work hst cross-references, and so forth. The 

information within database server 28 serves to facilitate storage and association of 
the image data files with one another, and to allow requesting clients to rapidly and 
accurately access image data files stored witiiin the system. Similarly, server 18 is 
coupled to one or more archives 30, such as an optical storage system, which serve 

25 as repositories of large volumes of image data for backup and archiving purposes. 

Techniques for transferring image data between server 18, and any memory 
associated with server 18 forming a short term storage system, and archive 30, may 
follow any suitable data management scheme, such as to archive image data 
following review and dictation by a radiologist, or after a sufficient time has lapsed 

30 since the receipt or review of the image files. 
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In the illustrated embodiment, other components of the PACS system or 
institution may be integrated with the foregoing components to fiirther enhance the 
system functionality. For example, as illustrated in Fig. 1, a 
compression/decompression library 32 is coupled to interface 20 and serves to store 
5 compression routines, algorithms, look up tables, and so forth, for access by 

interface 20 (or other system components) upon execution of compression and 
decompression routines (i.e. to store various routines, software versions, code tables, 
and so forth). In practice, interface 20 may be part of library 32. Library 32 may 
also be coupled to other components of the system, such as client stations 22 or 

10 printer interface 24, serving similarly as a library or store for the compression and 

decompression routines and algorithms. Although illustrated as a separate 
component in Fig. 1, it should be understood that library 32 may be included in any 
suitable server or memory device, including within server 18. Moreover, code 
defining the compression and decompression processes described below may be 

15 loaded directly into interface 20 and/or library 32, or may be loaded or updated via 

network links, including wide area networks, open networks, and so forth. 

Additional systems may be linked to the PACS, such as directly to server 28, 
or through interfaces such as interface 19. In the embodiment illustrated in Fig. 1, a 

20 radiology department information system or RIS 34 is hnked to server 18 to 

facilitate exchanges of data, typically cross-referencing data within database server 
28, and a central or departmental information system or database. Similarly, a 
hospital information system or HIS 36 may be coupled to server 28 to similarly 
exchange database information, workflow information, and so forth. Where desired, 

25 such systems may be interfaced through data exchange software, or may be partially 

or fully integrated with the PACS system to provide access to data between the 
PACS database and radiology department or hospital databases, or to provide a 
single cross-referencing database. Similarly, external chents, as designated at 
reference numeral 38, may be interfaced with the PACS to enable images to be 

30 viewed at remote locations. Such external clients may employ decompression 

software, or may receive image files aheady decompressed by interface 20. Again, 



15-IS-5393 



links to such external clients may be made through any suitable connection, such as 
wide area networks, virtual private networks, and so forth. 

Fig. 2 illustrates in somewhat greater detail the type of cross-referencing data 
5 made available to clients 20, 22, 24, 30 through database server 28. The database 
entries, designated generally by reference numeral 40 in Fig. 2, will include cross- 
referenced information, including patient data 42, references to specific studies or 
examinations 43, references to specific procedures performed 44, references to 
anatomy imaged 45, and further references to specific image series 46 within the 

10 study or examination. As will be appreciated by those skilled in the art, such cross- 

referenced information may include further information regarding the time and date 
of the examination and series, the name of diagnosing, referring, and other 
physicians, the hospital or department where the images are created, and so forth. 
The database will also include address information identifying specific images, file 

15 names, and locations of the images as indicated at reference numeral 48. Where the 

PACS includes various associated memory devices or short term storage systems, 
these locations may be cross-referenced within the database and may be essentially 
hidden from the end user, the image files simply being accessed by the system for 
viewing firom the specific storage location based upon cross-referenced information 

20 in the database. 

As described more fully below, in accordance with certain aspects of the 
present technique, descriptive information is used to identify preferred or optimal 
compression routines used to compreiss image data. Such descriptive information is 

25 typically available from header sections of an image data string, also as described in 
detail below. However, information available fi-om database server 28 may also 
serve as the basis for certain of the selections of the algorithms employed in the 
compression technique. Specifically database references may be relied upon for 
identifying such descriptive information as the procedures performed in an imaging 

30 sequence, specific anatomies or other features viewable in reconstructed images 

based upon the data, and so forth. Such information may also be available from the 
RIS 34 and fi-om tiie HIS 36. 
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Fig. 2 also illustrates an exemplary image file cross-referenced by the 
database entries. As shown in Fig. 2, image file 50 includes a plurality of image 
data sets 52, 54 and 56. In a typical image file, a large number of such image sets 
5 may be defined by a continuous data stream. Each data set may be compressed in 

accordance with specific compression algorithms, including lossless compression 
algorithms as described below, lossy compression algorithms, wavelet algorithms, 
and the preferred compression code table-based optimal compression algorithm 
described below. 

10 

Within each image data set, a descriptive header 58 is provided, along with a 
compression header 60. The headers 58 and 60 are followed by compressed image 
data 62. The descriptive header 58 of each data set preferably includes industry- 
standard or recognizable descriptive information, such as DICOM compUant 

15 descriptive data. As will be appreciated by those skilled in the art, such descriptive 

information will typically include an identification of the patient, image, date of the 
study or series, modahty of the system creating the image data, as well as additional 
information regarding specific anatomies or features visible in the reconstructed 
images. As described more fiilly below, such descriptive header data is preferably 

20 employed in the present technique for identification of optimal compression 

algorithms or routines used to compress the data within the compressed image data 
section 62. Data referring to the specific algorithm or routine used to compress the 
image data is then stored within compression header 60 for later reference in 
decompressing the image data. As described below, additional data is stored within 

25 tiie compressed image data, cross-referencing the algorithms identified in 
compression header 60 for use in decompressing the image data. Specifically, in a 
presently preferred embodiment, the compression header 60 includes identification 
of the length of subregions of the compressed image data, as well as references to 
specific optimal algorithms, in the form of compression code tables used to 

30 compress the subregions optimally. 
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Fig. 3 illustrates a typical image which is encoded by packets of digitized 
data assembled in a continuous data stream which may be compressed and 
decompressed in the present techniques. The image, designated generally by the 
reference numeral 100, will typically include features of interest 102, such as 
5 specific anatomical features. In medical diagnostic applications, such features may 

include specific anatomies or regions of a patient viewable by virtue of the physics 
of the image acquisition modality, such as soft tissue in MRI system images, bone in 
x-ray images, and so forth. Each image is comprised of a matrix having a width 104 
and a height 106 defined by the number and distribution of individual pixels 108. 

10 The pixels of the image matrix are arranged in rows 110 and columns 112, and will 

have varying characteristics which, when viewed in the reconstructed image, define 
the features of interest. In a typical medical diagnostic application, these 
characteristics will include gray level intensity or color. In the digitized data stream, 
each pixel is represented by binary code, with the binary code being appended to the 

15 descriptive header to aid in identification of the image and in its association with 

other images of a study. As noted above, such descriptive information may include 
industry standard information, such as DICOM comphant data. 

Fig. 4 graphically represents intensities of pixel data defining an image 
20 across a pair of rows of the image matrix. Each row of the image matrix will 

include a series of pixels, each pixel being encoded by binary data descriptive of the 
pixel characteristics, typically intensity. Thus, lighter regions of the reconstructed 
image will correspond to pixels having a higher intensity level, with darker regions 
having a lower intensity level. When illustrated graphically, the intensity levels of 
25 the pixels across the image may form a contour or trace as shown in Fig. 4. 

Specifically, Fig. 4 illustrates a first row 114 lying adjacent to a second row 116, 
each row including a series of pixels having various intensities as indicated by traces 
118 and 120, respectively. As will be appreciated by those skilled in the art, in 
practice, the graph of the pixel intensities would form a step-wise fimction along the 
30 position axis 122 and having amplitudes varying along the intensity axis 124. 
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It may be noted from Fig. 4 that in an actual image, variations in intensity 
along rows, and along columns as represented by correspondingly located intensities 
moving downwardly or upwardly through adjacent rows, will vary in accordance 
with the features represented in the reconstructed image. As illustrated in Fig. 4, 
5 row 1 14 includes areas of rising and declining intensity, including areas of both low 

intensity and high intensity. Row 116 includes areas of similar intensity, but 
varying by virtue of the features represented in the image. In accordance with the 
present technique, the image data stream is reduced to subregions as represented 
generally by reference numeral 126 in Fig. 4. While the subregions may be of 

10 different lengths (i.e. numbers of pixels), in the presently preferred embodiment, 

each subregion includes data encoding an equal number of pixels. Those skilled in 
the art will readily recognize, however, that after compression the actual length of 
codes for the subregion will vary depending upon the intensity of the pixels within 
the subregion and the dynamic range of the digital data encoding the pixel 

15 intensities. It should also be noted that where the row length of the image matrix is 

an integer multiple of the subregion width, individual subregions will align with one 
another moving down the image matrix as represented in Fig. 4. In general, 
however, the present technique is not limited to such integer multiple row widths. 

20 Each subregion of the image data stream may be analyzed to identify 

entropy levels of the image data for compression purposes. In general, the entropy 
of the image data refers to the relative variation in the pixel intensities within each 
subregion. Thus, although a specific subregion may include pixels of a high 
intensity or a low intensity, where the intensity level is relatively stable or constant 

25 in the subregion, entropy is considered to be low. Such regions are shown in Fig. 4, 

for example, in portions of traces 1 14 and 116 indicated by reference numerals 128 
and 132, respectively. On the contrary, where substantial variations in pixel 
intensity are encoded by the subregion data, the entropy level is considered high. 
Such regions are shown in Fig. 4 at reference numerals 130 and 134. It has been 

30 found that by subdividing the image data stream into subregions, which may have 

various lengths depending upon the image characteristics, and by analyzing the 
relative entropy level of each subregion, the subregions may be optimally 
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compressed in accordance with one of several candidate compression algorithms. 
Specifically, the present technique employs a series of predefined compression code 
tables which serve to translate pixel prediction errors into compressed codes for each 
subregion. The specific compression code table selected for each subregion is a 
5 function of the entropy identified for the subregion in the compression routine. 

As will be described more fiilly below, the present technique preferably 
employs pixel value predictors and identifies differences between predicted values 
of individual pixels (i.e. binary code for the pixel intensity or other characteristic) 

10 and the actual value for the respective pixels. In fact, several predictor algorithms 
may be employed, with the specific predictor being identified based upon image 
characteristics, such as characteristics encoded in a descriptive header for the image 
data stream. The predictors are based upon comparisons of target pixels, or pixels of 
interest, with neighboring pixels. Fig. 5 represents a pixel neighborhood which 

15 serves as the basis for references made in the predictor algorithms. The pixel 

neighborhood, identified by reference numeral 138 in Fig. 5, includes the pixel of 
interest 136, designated p(i, j). Neighboring pixels include a "northwest" pixel 140, 
designated p(i-l, j-1), a "north" pixel 142 designated p(i, j-1), a "northeast" pixel 
144 designated p(i+l, j-1), a "west" pixel 146 designated p(i-l, j). Some or all of 

20 these pixels may be employed in the predictor algorithm, or other additional pixels 

may be provided in a larger neighborhood. 

The preferred technique for compressing the image data stream in the system 
described above to create hybrid compressed image data files is summarized in Fig. 

25 6. The control logic may be subdivided into a series of logical blocks or segments, 

including a configuration segment 250 in which optimal compression algorithms, 
tables, predictor preferences, block or subregion length preferences, and so forth, are 
established. In a data access and retrieval segment 252 the image data to be 
compressed is received, as is descriptive information fi-om database server 28, if 

30 used. In an evaluation segment 254 the images are evaluated based upon the 

algorithms and preferences estabhshed in configuration segment 250. In a 
compression segment 256 certain values are computed and selected for the 
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compression routine, and the image data is compressed to create the resulting 
compressed file. File code is added to the compressed data in a final segment 258. 
The file is then completed at step 260 and stored at step 262. The compression and 
decompression logic is then completed by eventual access of the compressed image 
5 file at step 264, retrieval of a decompression algorithm or algorithms at step 266 and 

decompression at step 268 for viewing and output. The foregoing logical segments 
and processes will be described below in greater detail. 

The configuration segment 250 of the control logic includes steps of 
10 configuring the specific compression, predictor and block or subregion algorithms 

employed in the routine. Thus, at step 270 a series of compression code tables are 
generated for optimally compressing subregions of image data based upon relative 
entropy levels as indicated by prediction errors. As will be appreciated by those 
skilled in the art, such compression code tables serve to cross-reference original 
15 values with compressed values, generally in accordance with anticipated 

distributions or fi-equencies of occurrence of the original values. In the presently 
preferred embodiment, a series of compression code tables are estabUshed based 
upon analysis of typical images to be compressed by the routine. While the specific 
coding and ranges implemented in the compression code tables vary, and will 
20 generally be determined in specific apphcations empirically, examples of several 

such compression code tables are illustrated in Figs. 7-12. 

Referring to Figs. 7 - 12, each compression code table, such as table 170 in 
Fig. 7, comprises a series of compressed data values 172 cross-referenced to original 

25 image parameters 174. In the example illustrated in the Figures, analysis may 

include an indication of the resulting code links in bits, as indicated by reference 
nvimeral 176, and the span of code translated by each table entty, as indicated by 
reference numeral 178. In the illustrated embodiment, the compressed data code of 
column 172 translates difference values of column 174 as identified by application 

30 of predictors in the compression routine, followed by determination of the 

differences between the predicted and actual values. Thus, by way of example, 
where a difference value of 1 is identified based upon a predictor algorithm as 
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described below, table 170 provides a compressed data code of 100. Again, it 
should be noted that the difference value does not generally represent the value of 
the encoded pixel itself, but a difference between a predicted value and an actual 
value. Thus, in a code length of 3 bits, summarized in column 176 of table 170, the 
pixel of interest differing from the predicted value by 1 will obtain a compressed 
data value of 100, although the pixel value itself may be considerably longer. As 
also may be seen in Fig. 7, the compressed data codes provided in column 172 may 
be summarized as including a first portion 180 designatmg the level within the table, 
followed by a second portion 182 which indicates the position within the level or 
range. Thus, for example, in the case of the table of Fig. 7, a difference range of 
from -2 to -3 would be encoded as 1 101 followed by an additional bit set to a value 
of 0 or 1 depending upon whether the difference is -2 or -3. At an upper limit of the 
range, the compressed data code is taken to be the actual value of the individual 
pixel as described below with reference to Fig. 13. 

As can be seen from the Figures, the compression code tables for translating 
prediction errors or differences to compressed code are established to provide a 
range of coding appropriate to the levels or variations in the difference values for 
each subregion of the data stream. Specifically, table 170 of Fig. 7 is adapted to a 
lowest entropy level as indicated by the low difference variation (zero) 
accommodated by the shortest code in the table and the relatively fine steps between 
the difference levels. Fig. 8 represents a second compression code table 184 
providing for relatively higher entropy as indicated by the relatively broader 
maximum levels accommodated by the table and the relatively higher difference 
ranges as compared to table 170 of Fig. 7. Figs. 9, 10, 1 1 and 12 provide additional 
examples of compression code tables 186, 188, 190 and 192, respectively, for 
encoding successively higher entropy levels as indicated by prediction errors or 
differences within the various subregions. In the present embodiment illustrated, the 
code tables are constructed using a Huffinan code-based prefix and a multibit 
extension. 
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It should be noted that the compression ratio provided by each of the 
compression code tables in the family of tables varies depending upon the relative 
entropy level to be accommodated. Thus, table 170 of Fig. 7 provides a relatively 
high maximiun compression ratio of 16:1, while table 192 of Fig. 12 provides a 
lower maximum compression ratio of 2.29:1. However, it has been found that the 
use of a plurality of different compression code tables to accommodate the different 
entropies of the data stream subregions results in excellent compression due to the 
use of compressed data code lengths optimally selected to accommodate the 
variations in prediction differences for the subregions. 

It has been found that for specific types of images or for images having 
specific typical characteristics, various entropy levels may be anticipated. For 
example, in medical diagnostic imaging, relatively high entropy levels may be 
expected for specific modalities, such as CT and MRI data. Other imaging 
modalities may provide images with relatively lower variations in the image 
intensity levels, reflected by lower entropy values and correspondingly lower 
prediction differences. Moreover, it has been found that specific image types may 
provide higher or lower characteristic entropy values. In the medical diagnostics 
field, such image types may include specific anatomies, such as chest, head, 
extremities, and so forth, providing more or less variation, and stronger or weaker 
edge lines and contrast. The specific family of compression code tables, then, are 
preferably established based upon the typical images to be compressed by the 
system. 

Returning to Fig. 6, with the family of compression code tables generated 
and stored at step 270, configuration segment 250 continues with generation of 
predictor preferences at indicated at step 272. As described above, rather than 
employ actual pixel intensity values for encoding in the present compression 
technique, difference values may be employed based upon the predictor algorithms. 
The predictor algorithms generally produce relatively low values in low entropy 
regions of the image, and relatively higher values in higher entropy regions. 
However, the use of prediction errors for the data compression will generally result 
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in compression of values which are lower (i.e. shorter in length) than the original 
data stream code values due to the relatively gradual progression of intensity 
variations over the image. 

In the presently preferred embodiment, one or more of several predictor 
algorithms may be employed. Referring again to Fig. 5, in a simple and 
straightforward algorithm, the value of each pixel p(i, j) is predicted to be the value 
of the immediately preceding pixel p(i-l, j). This predictor algorithm provides a 
computationally extremely efficient tool for the prediction of each pixel value, with 
the first pixel in the image being estimated at a value of zero. The difference value 
is then generated by finding the absolute value of the difference between the 
predicted value and the actual value for the pixel of interest. The resulting 
difference values form a matrix of equal size to the original image matrix. 

Several ahemative predictor algorithms are presently contemplated, and 
others may be employed as well. In the presently preferred embodiment, these 
predictor algorithms may be summarized as follows: 



Ip(i, j) = ((3*ap(i-l J)) + Ip(i-1 J-1) + Ip(i, j-1) + Ip(i+1, j-l))/6 (P3); 



Ip(i,j) = Ip(i-l,j) 



(PI); 



Ip(i, j) = Ip(i-1, j) + Ip(i, j-1) - Ip(i-1, j-1) 



(P2); 



Ip(i, j) = ((-2*I(i-l, j-1)) + (3* Ip(i, j-1)) + (3*Ip(i-l, j)))/4 



(P4); 



Ip(i, j) = ((-5*I(i-l, j-1)) + (7* Ip(i, j-1)) + 
Ip(i+l,j-l) + (9*Ip(i-l,j)))/12 



(P5); 



Ip(i,j) = Ip(i,j-l) + Ip(i-l,j))/2 



(P6); 
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where the notation "I" represents the pixel intensity, and the individual pixel 
designations are mapped in accordance with Fig. 5. 

Again, it has been found that various predictors are particularly useful for 
5 various types of images, images resulting from various modalities, and various 

features or anatomies visible in the reconstructed image. For example, the predictor 
algorithm PI provides an extremely simple and rapid "last value" estimate of each 
current pixel value. Predictor P2 provides some degree of edge detection capability, 
but may accentuate the influence of noise in the image. The third predictor value P3 
10 is a linear predictor shown useftil on certain image types and employs values for 

more of the surrounding pixels, weighting the last value more heavily. In practice, 
the predictor preferences generated at step 272 of Fig. 6 will follow empirical results 
for various image types and modaUties. 

15 The final step illustrated in segment 250 of the control logic of Fig. 6 is the 

generation of subregion preferences for the division of the image data stream into 
subregions. As indicated above, while the lengths of subregions of the image data 
stream may be different or may vary depending upon such factors as the entropy 
level, in the presently preferred configuration, the subregion lengths are set equal to 

20 one another, but may vary depending upon such factors as the modality originating 

the image data, the image type, anatomy illustrated in the reconstructed image, and 
so forth. For example, subregion lengths of 32 pixels are presently set to a default, 
but alternative subregion lengths of 8, 16, 24, or more than 32 pixels, are 
contemplated. In general, the preferences set at step 274 in the present embodiment 

25 are dependent upon the number of columns in the image matrix, with subregion 
lengths of 8 being used for matrices having 64 columns or less, subregion lengths of 
16 being used for matrices having 128 columns or less, and so forth. To reduce the 
computational complexity and to improve speed, in the presently preferred 
embodiment, the subregions are taken along the major direction or axis (i.e. rows or 

30 colimms) of the image matrix, depending upon the direction in which the data is 

stored. Also, the subregions are preferably one dimensional only. 
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With the compression algorithms, tables, predictor preferences, and 
subregion size preferences set, the control logic illustrated in Fig. 6 advances 
through segment 252. In that segment of the logic, the image data is received as 
indicated at step 276. As noted above, the image data may be directly received from 
5 an extemal source, or may be accessed from a memory within the PACS itself At 

step 278, any descriptive information used to select algorithms or parameters of the 
compression routine are retrieved from the database server. As noted above, such 
descriptive information may be stored in databases, and may complement the 
information contained within the descriptive header of the image data stream. 

10 

Following segment 252, confrol advances to evaluation segment 254. 
Within this segment, the image data is reviewed for descriptive information as 
indicated step 280. As described above, where descriptive information is available, 
such as DICOM compliant data in a descriptive header section of the image data 

15 stream or descriptive data from the database, some or all of this data is reviewed at 

step 280. Based upon the preferences set in the configuration segment 272, 
predictors are selected at step 282 depending upon the image characteristics 
identified as step 280. Again, these may include the modality of the originating 
imaging system, the study type or anatomy featured in the image, the number of 

20 columns in the image, the number of rows, and so forth. Moreover, other factors 

may be considered in selecting the predictors at step 282, such as the computational 
efficiency desired, the processing power of the system, and so forth, with 
computationally efficient predictors being selected where such processor capabilities 
are limited, or where additional speed is desired. At step 284, the subregion size for 

25 division of the image data sfream into subregions is selected in accordance with the 
preferences estabhshed at step 274. Again, step 284 may consist of a default 
selection, which may be altered depending upon some or all of the characteristics or 
factors considered for selection of the predictors. 

30 Evaluation segment 254 continues with the selection of a subset of 

compression tables, where appropriate, as indicated at step 286. In particular, based 
upon certain image characteristics, it may be useful to preselect certain compression 
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tables as defaults. For example, specific image types originating in specific 
modalities, such as CT or MR images, may be best compressed with specific 
candidate tables which may be selected at step 286. At step 288 a compression 
header is affixed to the image data. As described below, this compression header 
5 will contain code identifying the version of the compression routine, the predictors 

selected at step 282, the subregion sizes selected at step 284, and so forth. 

At compression segment 256 of the control logic, a set of prediction errors or 
difference values are computed at step 290. As noted above, these values are based 
1 0 upon application of one or more of the predictor algorithms selected at step 282, and 

by subsequently comparing the predicted values with the actual values for each pixel 
to determine the prediction error or difference. The resulting prediction errors then 
form a data stream with the first pixel being assigned its actual value, followed by 
the difference values for each pixel of the subregion. 

15 

The present technique provides for optimal compression of each of the 
subregions based upon appropriate selection of tables fi-om the family of 
compression code tables established at step 270. To determine which tables provide 
best compression, at step 290 of the control logic of Fig. 6, candidate lengths of 

20 compressed data are computed by application of each of the candidate compression 

code tables to the data resulting firom application of the predictor algorithms. The 
total data stream length for each subregion is summed for each case (i.e. application 
of each of the compression code tables). Following completion of each subregion, 
the resulting sums are compared to identify which compression code table results in 

25 the shortest data stream for the subregion. The corresponding compression code 
table is then selected as indicated at step 292. At step 294, code identifying the 
selected table for each subregion is then inserted into the compressed data and its 
header, as described in greater detail below. 

30 As noted above, evaluation segment 254 or compression segment 256 may 

include defaults and constraints on the selection of the optimal code tables. For 
example, depending upon such factors as the bit depth of the image data to be 
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encoded, certain default selections among the compression code tables may be 
made, such as selection of the first four tables for bit depths of 8 or less. In addition, 
certain preferences for designated compression code tables may be made in this 
selection process, depending upon such factors as the modality originating the image 
5 data. By way of example, CT and MRI data may tend to encode higher transitions 

of pixel intensities, corresponding to higher entropy regions. Accordingly, 
preferences may be included in the selection of the compression code tables, 
depending upon the descriptive data reviewed at step 280, such as to prefer one or 
more higher entropy tables for images origiaating in such modalities. Finally, in the 

10 presently preferred embodiment, code will be inserted directly into the compressed 

data stream to provide an indication of which of the selected tables is employed for 
compressing the individual subregions. To optimize the compression, it may be 
desirable to limit the number of tables which may be selected in the preprocess 
segment of the control logic to a number which may be conveniently coded in a 

15 limited number of bits, such as 2 bits of data. Thus, where table identification codes 

are to be provided in 2 bits, a constraint may be imposed in the preprocess segment 
such that only four of the candidate tables may be selected for the compression. 
Where an additional bit is available for identification of the tables, 8 candidate tables 
may be selected and encoded. 

20 

At step 294 key code for the subregion compression is inserted mto the 
compressed data stream, immediately preceding the compressed subregion data. As 
noted above, where four candidate tables are employed for the compression, the 
code inserted at step 294 may include a pair of designated bits. The compression 
25 header created at step 288 cross references this key code to the selected compression 

tables. At step 296, the image data for the subregion is compressed by apphcation of 
the selected compression code table. The series of steps of segment 256 is repeated 
for each subregion of the image until the entire image is compressed. 

30 In final segment 258 the compressed data for each image is completed. 

Specifically, at step 298, a compression end block is added to the compressed data. 
At step 300, padding bits are inserted following the compression end block. Finally, 
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at step 302 a checksum value computed through the compression process is added to 
the compressed data to provide a means for verifying proper decompression. The 
position and type of code added to the compressed data during segment 258 is 
described more fully below with reference to Fig. 13. 

5 

Where a descriptive header is provided for the original image data, the 
descriptive header is preferably replaced adjacent to the compression header to 
complete the image file as indicated at step 260. It will be noted that the resulting 
data file is a hybrid compressed data file in which the descriptive header data is 

10 readable for image management, access, transmission, and similar purposes, while 

the bulk of the image data is optimally compressed as described above. At step 262, 
the resulting file is saved, typically in a short term storage system, or in an archive. 
At step 264 the image is accessed for reviewing, analysis, hard copy creation, and so 
forth. Upon access of the image, the decompression algorithms identified in the 

15 hybrid compressed data file are retrieved for decompression, and the compressed 

image data is decompressed, as indicated at step 266. This decompression follows 
the compression header information, with the same compression code tables being 
employed to decompress the data as were employed to compress it for each 
subregion. This application of the compression code tables results in identification 

20 of the prediction errors or differences on which the compression was based, and the 

specific predictor algorithm or algorithms employed to generate the differences are 
then used to regenerate the original image data in a lossless fashion. 

Fig. 13 represents an image data set compressed in accordance with the 
25 foregoing technique. The image data set, designated generally by the reference 

numeral 194, includes the descriptive header 196 appended at step 160, along with 
the compression header 198 affixed at step 192. Following the compression header 
198 is the compressed image data 200 which includes compressed image subregions 
202, 204 and so forth. Each compressed data subregion, in turn, includes an 
30 identification of the algorithm or algorithms (e.g. compression code tables) used to 

encode the subregion, followed by the actual compressed data. 
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Fig. 13 also illustrates a presently contemplated format for a compression 
header 198. As shown in Fig. 13, the compression header mcludes an identification 
of the compression routine version 206, followed by an identification of the 
predictors used in the compression process 207, and identification of the subregion 
5 length at reference numeral 208. A compression algorithm identification segment 

210 is then inserted in the compression header to identify which of the compression 
algorithms (e.g. compression code tables) were selected for the compression, and the 
manner in which each table is encoded in the data stream for each subregion. In the 
example of Fig. 13, for example, tables 0, 3, 4 and 5 were selected, as indicated at 
1 0 reference numeral 212, with each table being identified by a 2 bit binary code to be 

found within the first 2 bits of the compressed data for each subregion. 

Fig. 13 also illustrates the format in the presently preferred embodiment for 
the compressed data code for each image subregion 202, 204, and so forth. In the 

1 5 illustrated embodiment, a first portion of each subregion includes an identification of 

the compression table used to compress the subregion data, as indicated at reference 
numeral 214 for subregion 202, and reference numeral 218 for subregion 204. This 
key code is followed by the compression codes as indicated at reference numerals 
216 and 220, respectively. Finally, Fig. 13 illustrates the code inserted at the end of 

20 each compressed image file, or portion of a file. Specifically, following the last 

subregion BN of the image, the end of block code portion 222 is affixed. As noted 
above, this end of block code signals the end of an image, or may be used to signal 
the end of a portion of an image, where the data compression routine is changed 
within a single image, such as due to large variations in the image data enti-opy. The 

25 padding code is inserted as indicated at reference numeral 224. This code may be of 

variable size as needed to complete the compressed data file on a whole word length. 
Finally, a 32 bit checksum portion 226 is added to complete the compressed data. 

While the invention may be susceptible to various modifications and 
30 alternative forms, specific embodiments have been shown by way of example in 

the drawings and have been described in detail herein. However, it should be 
understood that the invention is not intended to be limited to the particular forms 
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disclosed. Rather, the invention is to cover all modifications, equivalents, and 
alternatives falling within the spirit and scope of the invention as defined by the 
following appended claims. For example, the foregoitig technique may be used on 
both images composed of var)dng gray levels, as well as on color images. As will 
5 be appreciated by those skilled in the art, color images will typically consist of 

various color components, which produce the appearance of color variations due 
to their respective intensities. The foregoing technique may be used either with or 
without separating the color components firom one another, but is preferably 
applied by separating the color components and processing (i.e. compressing) the 

10 components in groups. Similarly, multifirame images may be accommodated in 

the present technique. As will be appreciated by those skilled in the art, such 
images typically comprise a number of separate images encoded by rows and 
columns, without separate descriptive headers (e.g. DICOM compMant headers) 
positioned between the separate images. In such cases, code identifying offsets in 

15 the compressed data corresponding to locations of different image frames is 

preferably inserted into the descriptive header of the hybrid compressed data file. 
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CLAIMS: 

1 . A method for compressing image data from an uncompressed image 
data stream, the method comprising the steps of: 

(a) compiling a plurahty of compression mapping tables for converting 
5 uncompressed data representative of individual picture elements to compressed data; 

(b) applying at least first and second compression mapping tables from 
the plurality of compression mapping tables to subregions of an uncompressed 
image data sfream to compress the subregions; and 

(c) appending data for the compressed subregions to form a compressed 
1 0 image data sfream. 

2. The method of claim 1, wherein the compression mapping tables are 
compression tables mapping a parameter representative of each picture element to a 
compressed data code. 

15 

3. The method of claim 2, wherein the parameter is a prediction error 
for each picture element. 

4. The method of claim 3, wherein the prediction errors are identified 
20 by application of a desired predictor algorithm to the uncompressed image data 

sfream. 

5. The method of claim 1, including the fixrther step of selecting the 
compression mapping tables apphed in step (b) from the plurality of compression 

25 mapping tables. 

6. The method of claim 5, wherein the compression mapping tables are 
selected based upon relative enfropy levels of each subregion. 
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7. The method of claim 6, wherein the relative entropy levels are 
determined by analysis of relative variation of picture element intensity within each 
subregion. 

5 8. The method of claim 7, wherein the variation of picture element 

intensity is determined by application of a prediction algorithm to determine a 
difference between a predicted value of each picture element and the actual value of 
the respective picture element. 

10 9. The method of claim 1, wherein the compression mapping tables 

employed in step (b) are selected based upon which compression mapping tables 
provide the shortest compressed data stream for each subregion. 

1 0. The method of claim 1 , wherein the number of compression mapping 
1 5 tables employed in step (b) may be encoded with at most two bits of data. 

11. The method of claim 1, including the step of encoding in the 
compressed image data stream identifiers representative of the compression 
mapping tables applied in step (b). 

20 

12. A method for compressing image data, the method comprising the 

step of: 

(a) defining a family of compression code tables for converting 
uncompressed image data to compressed data; 
25 (b) storing the compression code tables in an image data compression 

station and in an image data decompression station; 

(c) selecting at least two of the compression code tables for compression 
of subregions of an image data stream; 

(d) compressing the image data stream in accordance with the selected 
30 compression code tables at the compression station for decompression at the 

decompression station. 
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13. The method of claim 12, including the step of encoding in the 
compressed image data stream identifiers of the selected compression code tables. 

14. The method of claim 12, wherein the compression code tables are 
5 defined based upon analysis of typical images to be compressed at the compression 

station. 

15. The method of claim 12, comprising the further step of applying a 
prediction algorithm to portions of the data stream representative of individual 

10 picture elements of an image to determine difference values between predicted 

values and actual values for the picture elements, and wherein the compression code 
tables are apphed to encode the difference values. 

16. The method of claim 12, wherein the compression code tables 
15 selected at step (c) are selected based upon which tables of the family of tables 

provides the shortest stream of compressed data for each subregion. 

17. The method of claim 12, wherein the number of compression 
mapping tables employed in step (c) may be encoded with at most two bits of data. 

20 

18. A system for storing, transmitting and viewing images, the system 
comprising: 

a data compression station configured to store a plurality of compression 
code tables for conversion of image data to compressed image data and to execute a 

25 compression routine in which an image data stream is converted to a compressed file 

by dividing into subregions and each subregion compressing in accordance with a 
compression code table selected from the plurality of compression code tables based 
upon which compression code table provides optimal compression of the subregion; 
a data storage device for receiving and storing the compressed file; and 

30 an image decompression station configured to store the plurality of 

compression code tables, to access the compressed file from the data storage device, 
and to execute a decompression routine in which the compression code tables 
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applied to compress the image data stream are applied to decompress the 
compressed file to reconstruct the image data stream. 

19. The system of claim 18, fiirther comprising a compression library for 
5 storing at least a portion of the compression and decompression routines, and 

wherein the compression station and the decompression station can access the 
compression library for code used in the compression or decompression routines. 

20. The system of claim 18, wherein the compression routine includes 
10 analysis of the image data stream for data representative of a characteristic of an 

image encoded by the image data stream. 

21 . The system of claim 20, wherein the characteristic is an identification 
of an image acquisition system originating the image data stream. 

15 

22. The system of claim 20, wherein the characteristic is an identification 
of a feature represented in an image encoded by the image data stream. 

23. The system of claim 18, wherein the compression routine includes 
20 encoding of identifiers of the selected compression code tables within the 

compressed file, and wherein the decompression routine includes analysis of the 
identifiers for selection of the same compression code tables for decompression of 
the compressed file. 

25 24. A computer program for compressing image data, the program 

comprising: 

a machine readable medium; and 

configuration code stored on the machine readable medium, the 
configuration code including an algorithm for analyzing an image data stream, for 
30 compressing subregions of the image data stream by application of a plurality of 
compression code tables, and for compiling the compressed subregions into a 
compressed data file. 
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25. The computer program of claim 24, wherein a family of candidate 
compression code tables is stored on the machine readable mediimi. 

5 26. The computer program of claim 24, wherein the algorithm includes 

computation of compressed data lengths provided by application of a plurahty of 
candidate compression code tables for compression of each subregion, and selection 
of the compression code tables providing the shortest compressed data lengths for 
each subregion. 

0 

27. The computer program of claim 24, wherein the code is installed on 
the machine readable medium via a configurable network link. 
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ABSTRACT OF THE DISCLOSURE 

An image data compression and decompression technique applies one or 
more compression code tables to optimally compress an image data stream. The 
compression code tables are established in accordance with anticipated image 
5 characteristics, and to accommodate different levels of variation or entropy in the 

image data. The unage data may be divided into blocks or subregions for analysis 
of which of the candidate compression code tables provides the optimal 
compression of each subregion. The appropriate code table is selected for each 
subregion. The evaluation of the compression performance based upon 
10 application of each compression code table may include analysis of prediction 

differences or errors between predicted values for pixels of an image and the 
actual values for the pixels. 
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As the bslow-named inventors, we hereby declare that: 

Our residences, post office addresses and citizenship are as stated below next to our names. 

We believe we are the original, first and joint inventors of the subject matter which is claimed and 
for which a patent is sought on the invention entitled: IMAGE DATA COMPRESSION EMPLOYING 
MULTIPLE COMPRESSION CODE TABLES, the specification of whicti is attached hereto. 

We hereby state that we have reviewed and understand the contents of the above-identified 
specification, including the claim, as amended by any amendment referred to above. 

We acknowledge the duty to disclose information which is material to the examination of this 
application in accordance with Title 37. Code of Federal Regulations, §1.S6(a). 



We hereby claim foreign priority benefits under Title 35, U.S.C. §1 19 of any foreign application for 
patent or inventor's certificate listed below, and have also identified below any foreign application for 
patent or inventor's certificate having a filing date before that of the application on which priority is claimed: 
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DATE OF FILING 
(day, month, year) 


PRIORITY CLAllWlED 
UNDER 37 U.S.C. 119 
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We hereby claim the benefit under Title 35, United States Code, §120 of any United States 
application(s) listed below and, insofar as the subject matter of each of the claims of this application is not 
disclosed in the prior United States application in the manner provided by the first paragraph of Title 35, 
United States Code, §112. 1 acknowledge the duty to disclose material information as defined in Title 37, 
Code of Federal Regulations, §1 .56(a) which occurred between ttie filing date of the prior application and 
the national or PCT international filing date of this application: 
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FILING DATE 
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both of GE IVledlcal Systems, 3000 North Qrandvlew Blvd., Waukesha, Wisconsin 5318B; Ronald E. 
Myrick (Reg. No. 26,315) and Henry J. Policinski, (Reg. 26,621) both of General Electric Company (W3E), 
3135 Easton Turnpike, Fairfield, Connecticut 06431-0001; and Patrick S. Yoder (Reg. No. 37,479). 
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patent and to transact all business in the Patent and Trademark Office connected therewith. 



Page 1 of 2 



NOV, 23, 1999 6:26PM GEMS IIS 847 704 8550 



NO. 7259 P. 3 



We hereby direct that all correspondence and telephone calls in connection with this application 
be addressed to Patrick S. Yoder. 7915 FM 1960 West, Suite 330, Houston. Texas 77070, (281) 970- 
4545, 

We hereby declara that all statements made herein of our own itnowledge are true and that all 
statements made on Information and belief are believed to be true; and further statements were made with 
the knowledge that willful false statements and the like so made are punishable by fine or imprisonment, 
or both, under Section 1001 of Title 18 of the United States Code and that all such willfiji fialse statements 
may jeopardize the validity of the application or any patent issued thereon. 
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